Evolutionary dynamics on degree-heterogeneous graphs 
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The evolution of two species with different fitness is investigated on degree-heterogeneous graphs. The 
population evolves either by one individual dying and being replaced by the offspring of a random neighbor 
(voter model (VM) dynamics) or by an individual giving birth to an offspring that takes over a random neighbor 
node (invasion process (IP) dynamics). The fixation probability for one species to take over a population of N 
individuals depends crucially on the dynamics and on the local environment. Starting with a single fitter mutant 
at a node of degree k, the fixation probability is proportional to k for VM dynamics and to l/k for IP dynamics. 
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In this letter, we investigate the likelihood for fitter mutants 
to overspread an otherwise uniform population on heteroge- 
neous graphs by evolutionary dynamics. Such a process un- 
derlies epidemic propagation 1 1], emergence of fads |2|], so- 
cial cooperation tj, or invasion of an ecological niche by a 
new species 0>BTIal- At each update event, two individuals 
from a total population N are chosen at random. One indi- 
vidual replicates while the other dies and is replaced by the 
newly-born offspring, so that N remains constant. A selective 
advantage, or fitness, also exists in which each each individ- 
ual may be a unit-fitness genotype 1 or genotype with lower 
fitness 1-5, with < s < 1. These fitnesses determine the 
replication or death rates of each individual. This selective ad- 
vantage leads to a dynamical competition in which selection 
dominates for large populations, while random genetic drift 
OlH occurs for small populations or weak selection. 

We consider three evolutionary models, distinguished by 
the order in which a pair individuals replicate and die: 
Biased Link Dynamics (LD) : A link is selected at ran- 
dom. If the individuals at the link ends are different, one 
of them is designated as the "donor" with probability pro- 
portional to its fitness. The replicate of the donor then re- 
places the other individual: 10 — » 11 with probability 1/2 
while 10 -> 00 with probability (l-s)/2 (Fig.[l). 
Biased Voter Model (VM) : An individual dies with prob- 
ability inversely proportional to its fitness, and is then replaced 
by the offspring of a randomly-chosen neighbor. Equivalently, 
death occurs randomly and replacement is proportional to the 
fitness of the donor. We implement the VM by updating a 
randomly-chosen genotype with probability 1, while the fit- 
ter genotype 1 is updated with a probability 1-5. Each indi- 
vidual in this death-first/birth-second process can equivalently 
be viewed as a voter that adopts the opinion of a randomly- 
selected neighbor mini in. 

Biased Invasion Process (IP) : In this birth- 
first/death-second process, a randomly-chosen individual 
replicates with probability proportional to its fitness, and its 
offspring then replaces an individual at a randomly-chosen 
neighboring node 0|5l01- 

One genotype ultimately replacing all other genotypes in 
the population is termed fixation. An important, and easily 
checked fact is that these evolutionary models are equivalent 
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FIG. 1: Update illustration for two specific nodes. Genotypes and 
1 are denoted by O and • respectively. Shown are the possible tran- 
sitions and their respective relative rates due to the interaction of two 
nodes across a link for LD, VM, and IP dynamics. 

on degree-regular graphs; moreover, as we will show, the fix- 
ation probability in LD can be obtained exactly, independent 
of the underlying graph. However, essential differences arise 
on degree-heterogeneous networks flll[l2l[l3ll that may lead 
to an enhancement of the fixation probability, as discovered 
previously for the IP |@|. Here we cast LD, VM, and IP on 
degree-heterogeneous graphs within the same unifying frame- 
work to understand the interplay between selection and ran- 
dom drift on the fixation probability. By this approach, we 
show that on degree-heterogeneous graphs the best strategy to 
reach fixation with VM dynamics is for the fitter genotype to 
be on high-degree nodes. Conversely, for IP dynamics, it is 
best for the fitter genotype to be on low-degree nodes. 

We first study the evolution in VM dynamics. We symbol- 
ically represent the state of the system by r\. In an elemental 
time interval St we choose a random node x. If the geno- 
type at this node at time t, denoted as rf{x), equals 0, then 
node x is updated by choosing a random neighbor y and set- 
ting rf t+5, \x) = rf(y) (Fig. [I). However if 77' (x) = 1, the VM 
update is implemented with probability 1 - s. This update rule 
can be written as 

Pfo-^*] = J] ^-r l^-n(x)]T,(y) + (l-s)T}(x)[l-T,(y)]) (1) 

where r/ x denotes the state obtained from 77 by changing only 
the genotype at node x. The first term describes the update 
step for the case where (r/(x), rj(y)) = (0, 1) and x,y are con- 
nected. Each of the nearest neighbors y of x may be selected 
with probability A xy /k x . Here A xy is the adjacency matrix 
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whose elements equal 1 if xy are connected and zero other- 
wise. The second term in Eq. Q is explained analogously. 

For degree-heterogeneous graphs, the density p k of geno- 
type 1 at nodes of degree k increases by 1 /N k with probability 
Fk{r]) and decreases by 1 /N k with probability B k (rj) in an ele- 
mental update, where 
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are the forward (0 — > 1) and backward (1 — » 0) evolution rates. 
The primes on the sums denote the restriction that the degree 
of nodes x equals k. The sum over all k then gives the total 
transition rate of Eq. Q. We seek the fixation probability O 
to the state consisting entirely of genotype 1 as a function of 
the initial densities of 1. This probability obeys the backward 
Kolmogorov equation GO = PhIi . subject to the boundary 
conditions O(0) = and O(l) = 1. In the diffusion approxi- 
mation, the generator G of this equation may be expressed as 
a sum of the changes in p k over all k, 
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with 5p k = 1 /N k = 1 /(N n k ) the change in p k in a single update 
of a node of degree k, and d k = S-. 

b ' " dpi 

For the special case of degree-regular graphs, where k x — k 
for all nodes, both sums in Eq. count the total number a of 
active links between different genotypes 
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where the moments of the degree distribution are defined by 
p„ = N~ l Yjx k" - Yik k"n k . The generator thus reduces to 



G = a 



(5) 



where we use dp = 6t = 1 /N. In this form, the convection 
and diffusion terms differ by a factor 0{sN). Thus selection 
dominates when the population N is larger than 0(1/5), while 
random genetic drift is important otherwise. 

Notice that the probability of increasing the density of 
genotype 1 at each update is a factor 1/(1 - s) larger than 
the probability of decreasing the density. By its construc- 
tion, this same bias arises for LD on general networks. As 
a consequence of this bias, the evolutionary process underly- 
ing fixation is the same as the absorption of a uniformly biased 
random walk in a finite interval 11411 . from which the fixation 
probability is 
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The former is the exact discrete solution of GO = on a finite 
network, while the latter continuum limit represents the solu- 
tion to GO = in the diffusion approximation. These results 
apply for all three models on degree-regular graphs and for 
LD on general graphs. 

For degree-heterogeneous graphs, the conserved quantity 
for neutral dynamics {s — 0) is the average degree-weighted 
density a>\ II llll 211 . where the degree-weighted moments are 
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while the overall density p of genotype 1 is no longer con- 
served. The existence of this new conservation law suggests 
that we study the time evolution of the expectation value of 
oj] which we henceforth denote as to for notational simplicity. 
Since 01(77.1) = to(r/) + k x (l - 2rj(x))/p]N, 
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Notice that a> is conserved in the absence of selection (s = 0) 
a feature that ultimately stems from the update rate being in- 
versely proportional to node degree [Eq. Q]. To evaluate 
the expression in Eq. © we make the mean-field assump- 
tion that the degrees of connected nodes in the graph are un- 
correlated. Thus we replace the elements of the adjacency 
matrix by their expected values, A„ = k x k y /p\N. This as- 
sumption simplifies Eq. (|8|i to d t u> = su>(l - of), with solution 

cj(ty l = 1 - [i - Gj(0)- l ] e - s ' . 

For uncorrelated graphs, Eqs. (|2} simplify to 
^kiv) = n k cj{\ - p k ), ~R k {tj) = (1 - s)n k {\ - oj)p k . (9) 
Thus the time evolution of the expectation value of p k is 
S Pk (F k -B k ) 



d,pk 
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To solve this equation we combine it with d,<±> = sw(l - ai) to 
give d,(a) — pk) = -(oj - p k )(l - 5(1 - to)), with solution 

cj(f)-p k (f) = e -'M0)-p A .(0)]{w(0) + [l-oj(0)]e- st ] . (11) 

For small selective advantage (s <k 1), this equation involves 
two distinct time scales. On a time scale of order one, all the 
p k become equal to a>, whereas the evolution of ai occurs on a 
longer time scale of order s~ l » 1 (Fig. |2}. 

We now determine the fixation probability simply by re- 
placing the p k by a> in the forward and backward rates F and 
B in Eqs. @. In a similar vein, we replace the derivative d k 
by (kn k jp.\)d ul ITlll . Then the generator in Eq. becomes 
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FIG. 2: Moments of the 1 density in the biased VM and biased IP on a 
network of 10 4 nodes with a power-law degree distribution n k ~ k~ v 
(v = 2.5), and no correlations between node degrees. Nodes with 
degree larger than the mean degree are initialized to 1 while all other 
nodes are 0. For the VM, s = 8.5 x 10~ 4 , while for the IP, s = 10~ 4 . 



FIG. 3: Scaling plot of fixation probabilities for VM (filled) and IP 
dynamics (open symbols). Data are for degree-uncorrelated graphs 
with N = 10 3 , /Ji = 8, and degree distribution exponent v = 2.1 (O), 
2.5 (A), or 3.0 (□). Initially each node is a mutant with probability 
1/2, (ll>i = LO-i = 1/2). The curve corresponds to Eqs. J 141 or <19> . 



Apart from an overall constant for the time scale, this gen- 
erator is identical to that of degree-regular graphs l|5} if we 
replace /V by an effective population size N s s = N^\ln2- For 
a network of /V nodes with a power-law degree distribution, 

n k oc k~ v , N e fi scales as fllll 
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with logarithmic corrections for v = 2 and v = 3. Thus N e s 
becomes much less than /V when \i% diverges; this occurs when 
v < 3. A similar change in the effective size of the population 
is observed for biological speci es evolving in a spatially het- 
erogeneous environment I4 ll5ll . 

The solution to GO = 0, with G given by Eq. H 2t is 
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Our numerical data for the fixation probability shows both ex- 
cellent scaling and agreement with this functional form for <D 
(Fig- 13- Eq. 1141 also provides the fixation probability when 
the system starts with a single mutant at a node of degree k: 
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The crucial feature is that the fixation probability of a single 
fitter mutant is proportional to the degree of the node that it 
initially occupies (Fig. [4}. Notice also that because the relative 
effect of selection versus random genetic drift is determined 
by the variable combination sN e s, random genetic drift can be 
important for much larger populations compared to the case 
of degree-regular graphs. In fact, for a power-law graph with 
v < 2, random genetic drift prevails for all population sizes. 

We now study fixation in the complementary biased inva- 
sion process. Here a randomly-selected individual reproduces 
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FIG. 4: Fixation probability of a single mutant initially at a node 
of degree k on an uncorrelated power-law degree distributed (n k ~ 
kT v , v = 2.5) graph with N = 10 3 and = 8. The empty symbols 
correspond to IP dynamics with s = 0.004 (□), s = 0.008 (O) and 
s = 0.016 (A); the filled symbols correspond to VM dynamics with 
s = 0.01, (■), s = 0.02 (•) and s = 0.08 (a). The solid lines, with 
slopes +1 and -1, correspond to the second of Eqs. I15i and I2UI . 



with probability proportional to its fitness; hence the transition 
probability is 

Pto^7J=£^{[i-^*)]'?(y)+(i-J)'KJf)[i-»7()')]).(i6) 

y 

Notice an essential difference between VM and IP dynamics. 
In the VM the transition rate is proportional to the inverse 
degree k x of the node of the disappearing genotype [Eq. flj], 
while in the IP the transition rate is proportional to the inverse 
degree k y of the node of the reproducing genotype [Eq. ( I16H . 

For degree-uncorrelated graphs, the transition probabilities 
are 

k k(l - s) 
Fjt(77) = —n k p(l -p k ) , Bjfc(?7) = n k (l -p)p k ■ (17) 

Consequently the time evolution of p k is given by, in analogy 
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with Eq. ( 1 10t . d,pk = -j^[p - pk + spuil - p)], from which 
low-order moments obey the equations of motion: 
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In contrast to the VM, the conserved quantity in the unbiased 
IP is co-i, the inverse degree-weighted frequency. For the bi- 
ased IP, co-i becomes the most slowly changing quantity (see 
Fig. |3 ■ Hence we transform all derivatives with respect to pk 
in the generator to derivatives with respect to <y_i to yield 
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from which, in close analogy with our previous analysis of the 
VM, the fixation probability is 
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From Eq. d!9i . the effective population size iV e fr equals N, 
contrary to VM dynamics (Eq. dl4» . More strikingly, the fix- 
ation probability of a single mutant acquires the non-trivial 
dependence of the degree k of the occupied node (Fig. [4} 
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To conclude, mutants are more likely to fixate in the voter 
model (VM) when they are initially on high-degree nodes 
[Eq. (I15H . while in the invasion process (IP) fixation is more 
probable when mutants start on low-degree nodes [Eq. (I20H . 
This behavior is understandable simply. In the VM, a well- 
connected individual is more likely to be asked his opinion 
before he asks one of his neighbors. In the IP, a mutant on a 
high-degree node is more likely to be invaded by a neighbor 
before the mutant itself can invade. Thus network heterogene- 
ity leads to effective evolutionary heterogeneity. 

We can also understand the evolution when a mutant ap- 
pears at a random node on a graph. In the selection-dominated 
regime (sN e s » 1) of the VM, we average Eq. dl5> over 
all nodes and find that the fixation probability on degree- 
uncorrelated graphs is smaller by a factor n\ln2 ^ 1 than that 



on regular graphs. Thus a heterogeneous graph is an inhos- 
pitable environment for a mutant that evolves by VM dynam- 
ics. Conversely, performing the same average of Eq. d20l over 
all nodes, the fixation probability for the IP is the same on 
all degree-uncorrelated graphs. Finally, in the small-selection 
limit (sN e « 1), the node average fixation probability is the 
same for both the VM and IP on degree-uncorrelated graphs. 
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